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Semileptonic decays of heavy baryons consisting of one heavy {Q = b, c) and two 
hght {q = n, d, s) quarks are considered in the heavy-quark-hght-diquark approxi- 
mation. The relativistic quasipotential equation is used for obtaining masses and 
wave functions of both diquarks and baryons within the constituent quark model. 



^^ , The weak transition matrix elements are expressed through the overlap integrals of 



the baryon wave functions. The Isgur-Wise functions are determined in the whole 

accessible kinematic range. The exclusive semileptonic decay rates and different 

""Oj' asymmetries are calculated with applying the heavy quark l/mq expansion. The 

J.. . evaluated A^ — > Ac/i/ decay rate agrees with its experimental value. 
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I. INTRODUCTION 

The description of heavy baryon properties represents a very interesting and important 
problem in quantum chromodynamics. Since the baryon is a three-body system, its theory 
is much more complicated compared to the two-body meson system. The quark-diquark 
picture of a baryon Jl|J^ is the popular approximation widely used to describe the baryon 
properties [l|, |2, |3|, y, |5|. Such approximation allows to reduce the very complicated rela- 
tivistic three-body problem to the two-body one. Recently, we evaluated the masses of the 
ground state heavy baryons in the framework of the relativistic quark model based on the 
quasipotential approach p. The heavy-quark-hght-diquark picture of the heavy baryons 
was assumed. Both scalar and axial vector light diquarks were considered. The relatively 
large size of the light diquark was effectively taken into account by calculating the diquark- 
gluon interaction form factor through the overlap integral of the diquark wave functions. All 
the parameters of the quark model had fixed values which were determined from the previous 
studies of heavy and light meson properties [7|, |8|, |9|, |lO| . The overall reasonable agreement 
(within a few MeV) of our model predictions for heavy baryon masses with the available 
experimental data supplies further support for the use of the heavy-quark-light-diquark 
approximation. 



In this paper we continue the study of heavy baryon properties and apply our relativistic 
quark model for the calculation of their exclusive semileptonic decays. Such investigations 
are important, since they provide an additional source (complimentary to the heavy me- 
son weak transitions) for determining the parameters of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix, such as Vcb, from the comparison of the theoretical predictions with the 
experimental data. We limit our present consideration to the heavy-to-heavy [b —>■ c) tran- 
sitions, where both the initial and final baryons contain heavy quarks. For such transitions 
the heavy quark effective theory (HQET), which is based on the l/mg expansion of the 



QCD Lagrangian and the emerging heavy quark symmetry [11|, provides the most effective 
constraints on theoretical models and significantly reduces the number of independent form 
factors in each order of the heavy quark expansion. In particular, transitions involving the 
^Q {Q = b, c) baryons have the simplest structure, since the spectator light degrees of free- 
dom (the scalar diquark) for these baryons have zero angular momentum. ^ In the heavy 
quark limit only one universal form factor, the so-called Isgur-Wise function, is required 
to describe the Aq — »• Aq' transition [ij, [l3[ . At the subleading order of the heavy quark 
expansion one mass parameter and one additional function emerge [ij]. The consideration 
of the Qq baryon decays is considerably more complicated, since light degrees of freedom 
(the axial vector diquark) now have spin 1. It is necessary to introduce two functions for 
parameterizing the VIq -^ Q,q; transition in the heavy quark limit 12|, and five additional 



functions and one mass parameter are needed at the subleading order in I/ttiq [15|- Note 
that transition matrix elements between baryons with spectator diquarks having different 
spins (e.g., Af, — »• S,,, which violate isospin symmetr y) v anish in the heavy quark limit and 
can proceed only due to the subleading corrections [12|- The heavy quark symmetry alone 
does not allow the determination of the corresponding Isgur-Wise functions and mass pa- 
rameters. Only the normalization of some of these functions is known at the point of zero 
recoil of the final heavy baryon. Thus, for the determination of these functions in the whole 
kinematic range the application of nonperturbative methods is necessary. Many different 
approaches were previously used for the calculation of the Isgur-Wise functions of heavy 
baryons. However, most of them have important limitations. In some of these approaches 
the Isgur-Wise functions are calculated only at one kinematic point or in a limited region 
and then extrapolated to the whole kinematic range using an ad hoc ansatz, while other ap- 
proaches assume some parameterization for the heavy baryon wave functions. The main aim 
of this paper is to determine the corresponding Isgur-Wise functions in the whole kinematic 
range through the overlap integrals of the heavy baryon wave functions in a consistent way 
within the relativistic quark model. These wave functions are known from the previous cal- 
culation of baryon masses p . On this basis exclusive semileptonic decay rates and different 
asymmetries can then be evaluated within the heavy quark expansion. 

The paper is organized as follows. In Sec. |n] we describe our relativistic quark model 
and present predictions for the masses of ground-state light diquarks and heavy baryons in 
the heavy- quark-light diquark picture. In Sec. Illll we discuss the determination of the weak 
current matrix element between heavy baryon states. The relativistic transformation of the 
baryon wave function from the rest to the moving reference frame is presented. The general 
expressions for the weak matrix elements, decay rates and different asymmetries for heavy 
baryons with scalar and axial vector diquarks are given in Sec. IIVI In Sec. |3 semileptonic 



^ The structure of the decay matrix elements for the Aq baryons is simpler than for heavy mesons, since in 
the latter case light degrees of freedom have spin 1/2. 



decays of heavy baryons with the scalar diquark are considered using the heavy quark 
expansion. Exphcit expressions for the leading and subleading Isgur-Wise functions are 
obtained as the overlap integrals of the baryon wave functions. The predictions for decay 
rates and the slope of the Isgur-Wise function are compared with the experimental data 
for the Afo —>■ AcCu decay. Semileptonic decay rates of heavy baryons with the axial vector 
diquark are studied in Sec. IVII in the heavy quark limit. Finally, the comparison of our 
results for the heavy baryon semileptonic decay rates with previous theoretical predictions 
and our conclusions are given in Sec. IVIII 



II. RELATIVISTIC QUARK MODEL FOR HEAVY BARYONS 



In the quasipotential approach and quark-diquark picture of heavy baryons p the inter- 
action of two light quarks in a diquark and the heavy quark interaction with a light diquark 
in a baryon are described by the diquark wave function {"^d) of the bound quark-quark state 
and by the baryon wave function (\I/b) of the bound quark-diquark state respectively, which 
satisfy the quasipotential equation of the Schrodinger type 

^'''^''' P'UMP)-f^v(P,r,,M)^Mn), (1) 
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where the relativistic reduced mass is 



_ E,E2 _ M^ - {ml - mlf 
^^ " El + ^2 " 4M3 ' ^^> 



and El, E2 are given by 



M^ — 1711 + mj 77. ^^ — m1 + ml 

~2M ' ^~ 2iir 
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Here M = Ei + E2 is the bound state mass (diquark or baryon), ?7ii 2 are the masses of 
light quarks (gi and ^2) which form the diquark or the masses of the light diquark (d) and 
heavy quark {Q) which form the heavy baryon (B), and p is their relative momentum. In 
the center of mass system the relative momentum squared on mass shell reads 

.2. ,.^ _ [M' - (^1 + m^Ym^ - (mi - m2)^] 

The kernel V{p, q; M) in Eq. (^ is the quasipotential operator of the quark-quark or 
quark-diquark interaction. It is constructed with the help of the off-mass-shell scattering 
amplitude, projected onto the positive energy states. In the following analysis we closely 
follow the similar construction of the quark-antiquark interaction in mesons which were ex- 
tensively studied in our relativistic quark model |7i]. For the quark-quark interaction in a 
diquark we use the relation Vqg = Vqq/2 arising under the assumption about the octet struc- 
ture of the interaction from the difference of the qq and qq colour states. An important role 
in this construction is played by the Lorentz-structure of the confining interaction. In our 
analysis of mesons while constructing the quasipotential of the quark-antiquark interaction, 
we adopted that the effective interaction is the sum of the usual one-gluon exchange term 
with the mixture of long-range vector and scalar linear confining potentials, where the vector 



confining potential contains the Pauli terms. We use the same conventions for the construc- 
tion of the quark-quark and quark-diquark interactions in the baryon. The quasipotential 
is then defined by [6,, J,, J^] 

(a) for the quark-quark (qq) interaction 



Vqq{p, q; M) = Ui{p)u2i-p)Vqq{p, q] M)ui{q)u2{-q), 



(5) 



V,,(p,q;M) = - 



^a,D,Ml'il2 + kL (k)r^(k)r2;,(-k) + KL(k) 



(b) for quark-diquark {Qd) interaction 

l^Q.(p,q;M) = M^E2MnQ(p)|a.D,.(k)r«Q(g) 
2^Ed{p)Ed{q) 3 

+ra{P)uQ{p)VgA^)uQ{q)UQ). 



(6) 



where as is the QCD coupling constant; {d{P)\J^\d(^)) is the vertex of the diquark-gluon 
interaction which takes into account the diquark size 'G^ in terms of the diquark wave function 
overlap 



gluon propagator in the Coulomb gauge 

A-n 



P = {Ea, -p) and Q = [E^, -q), E^ = (M^ - m^ + Mj)/(2M) J. D^,{k) is the 



I^-(k) = -^, D^^(k) = -^ 5--^ 



and k = p — q; 7^ and u{p) are the Dirac matrices and spinors 
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s(p) + m 



with e{p) = Vp^ + m'^. 

The diquark state in the confining part of the quark-diquark quasipotential (jHI) is de- 
scribed by the wave functions 
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^d{p) 



for scalar diquark 

for axial vector diquark 



where the four vector 



£dip) 



'i^dP) 



£d + 



(£dP)p 



(9) 



(10) 



is the polarization vector [s'^{p)p^ = 0] of the axial vector diquark with momentum p, 



Ed{p) = yp^ + Mj and £^(0) = (0, £d) is the polarization vector in the diquark rest frame. 
The effective long-range vector vertex of the diquark can be presented in the form 



J,- 



[P + Q), 



d]^ 



2JE,{p)E,{q) 



jP + Q), 



2JEd{p)E,{q) 



+ 2Md t"^ 



for scalar diquark 

for axial vector diquark 
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where k = (0,k). Here the antisymmetric tensor reads 

{J:,^)l = -z{g,,6:-g,j;), (12) 

and the spin S^ of the axial vector diquark is given by {Sd-k)ii = —i^ku- We choose the total 
chromomagnetic moment of the axial vector diquark fid = llll] • Such a choice appears to 
be natural, since the long-range chromomagnetic interaction of the diquark proportional to 
Hd then vanishes in accord with the flux tube model. 

The effective long-range vector vertex of the quark is defined by [7|, |]J 



T,{k)=j, + ^a,,~k'', fc = (0,k), (13) 

where k, is the Pauli interaction constant characterizing the anomalous chromomagnetic 
moment of quarks. In the configuration space the vector and scalar confining potentials in 
the nonrelativistic limit reduce to 

KonfW = (l-^)Konf(r), 

VLfir) = eKonf(r), (14) 

with 

Konf (r) = Konf (0 + Konf (0 = Ar + B , (15) 

where e is the mixing coefficient. 

The constituent quark masses rrib = 4.88 GeV, rric = 1.55 GeV, m„ = m^ = 0.33 GeV, 
rris = 0.5 GeV and the parameters of the linear potential A = 0.18 GeV^ and B = —0.3 
GeV have the usual values of quark models. The value of the mixing coefficient of vector 
and scalar confining potentials e = —1 has been determined from the consideration of 
charmonium radiative decays [8| and the heavy quark expansion [9|. Finally, the universal 
Pauli interaction constant k = —1 has been fixed from the analysis of the fine splitting of 
heavy quarkonia ^Pj- states j3]. Note that the long-range chromomagnetic contribution to 
the potential in our model is proportional to (1 + k) and thus vanishes for the chosen value 
of K = —1. 

The quasipotential (0) can be used for arbitrary quark masses. The substitution of the 
Dirac spinors into (0) results in an extremely nonlocal potential in the configuration space. 
Clearly, it is very hard to deal with such potentials without any additional transformations. 
In oder to simplify the relativistic qq potential, we make the following replacement in the 
Dirac spinors 



Illy tm 

,l3,Ii3 



ei,2(p) = JrnJ2 + P^ -> ^i,2- (16) 



This substitution makes the Fourier transformation of the potential (0) local, but the re- 
sulting relativistic potential becomes dependent on the diquark and baryon masses in a 
very complicated nonlinear way. We consider only the baryon ground states, which further 
simplifies our analysis, since all terms containing orbital momentum vanish. The detailed 
expressions for the relativistic quark potential can be found in Ref . p . The obtained masses 
of the light diquarks are given in Table IH The heavy baryon masses calculated in the heavy- 
quark-light diquark approximation are presented in Table |n]in comparison with the available 
experimental data [l9| . There an overall good agreement of our predictions with experiment 
is found. 



TABLE I: Masses of light ground state diquarks (in MeV). S and A denote scalar and axial vector 
diquarks which are antisymmetric [q, q'] and symmetric {q, q'} in flavour indices, respectively. 



Quark content 



Diquark type 



Mass 



[u,d] 
{u,d} 

[u,s] 
{u,s} 
{s,s} 



S 
A 
S 
A 
A 



710 
909 
948 
1069 
1203 



TABLE II: Masses of the ground state heavy baryons (in MeV). 



Baryon 


i{J' 


j^^theor ^ 


M^^P [19] 


Ac 


o(r; 


2297 


2284.9(6) 


Se 


i(r: 


2439 


2451.3(7) 


K 


i(r; 


2518 


2515.9(2.4) 


^c 


^(r 


) 2481 


2466.3(1.4) 


^c 


1(1+ 


) 2578 


2574.1(3.3) 




1(3 + 
2V2 


) 2654 


2647.4(2.0) 


Qc 


o(r; 


2698 


2697.5(2.6) 


ni 


o(r: 


2768 




Ab 


o(r; 


5622 


5624(9) 


Sft 


i(r: 


5805 




s* 


i(r; 


5834 




Sfe 


1(1+ 

2V2 


) 5812 




"b 


1(1+ 

2V2 


) 5937 




"* 
^b 


1(3 + 
2V2 


) 5963 




^b 


o(r: 


6065 




ni 


o(r; 


6088 





III. MATRIX ELEMENTS OF THE WEAK CURRENT FOR HEAVY BARYON 

DECAYS 



In order to calculate the exclusive semileptonic decay rate of the heavy baryon, it is 
necessary to determine the corresponding matrix element of the weak current between baryon 
states. In the quasipotential approach, the matrix element of the weak current J^ = 
(5'7/x(l — l5)Q, associated with the heavy-to-heavy quark Q ^ Q' [Q = b and Q' = c) 
transition, between baryon states with masses Mbq , Mb^, and momenta pbq , Pb„, takes the 
form 12 



iBQ'iPB^,)\j;:\BQiPB^)) 



(27r)6 



^Bo, PB,, (p)rM(p, q)*BQ PBo (q)' (i7) 



^Q' ^Bq^ 



where r^(p, q) is the two-particle vertex function and ^ b-pb ^i'^ ^^^ baryon {B = Bq, Bqi) 
wave functions projected onto the positive energy states of quarks and boosted to the moving 
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FIG. 1: Lowest order vertex function T^^' contributing to the current matrix element (|17j) . 
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FIG. 2: Vertex function F^^^ taking the quark interaction into account. Dashed hues correspond 
to the effective potential Vqd in ©• Bold lines denote the negative-energy part of the quark 
propagator. 



reference frame with momentum pe- 

The contributions to F come from Figs.^andEJ The contribution F*^^) is the consequence 
of the projection onto the positive-energy states. Note that the form of the relativistic 
corrections resulting from the vertex function F^^^ is explicitly dependent on the Lorentz 
structure of the quark-diquark interaction. In the heavy quark limit rriQ -^ cxd only F*^^) 
contributes, while F*^^^ gives the subleading order contributions. The vertex functions are 
given by 



ri'Hp,q) = V^d(Prf)^Q'(PQ')7^(l - l'')uQ{qQ)i;d{qd){2Tif6{^d - Qd) 



(18) 



and 



F(f)(p,q) = rAPd)uQipQ>){^,{l-l' 



{-) 



A^o'ik) 



-7°VQd(prf - Qd) 



+VQd(pd - qd) ^ ,(k')'+e i{ — \^^^^^^ " 7^)}^Q(gQ)^rf(gd) 



(19) 



where the superscripts "(1)" and "(2)" correspond to Figs. C] and |21 k = pq/ — A; k' 
qg + A; A = Mo v' - Mb.v; e(p) = Vm^ + p^; 



tS-\v) 



e{p) - (m7° + 7°(7p)) 



Here |2L. 

3 



PQ',d = eQ/,rf(p)t;'±^nW(t;>% v 






) 



^1 ^B,, 

3 /^ 

gQ,d = eQ,,(g)t; ± ^n«(t;)g\ t;^ = ^, (20) 

Mbo 



i=l 



and n*^*^ are three four-vectors given by 



s,J 



V V' 



It is important to note that the wave functions entering the weak current matrix element 
(fTTjl are not in the rest frame in general. For example, in the Bq baryon rest frame (v = 0), 
the final baryon is moving with the recoil momentum A. The wave function of the moving 
baryon ^E'b^, a is connected with the wave function in the rest frame "^Bq/ o = "^Bq, by the 



transformation |2_, 



^B,,^{p) = D'ci'{RfjDJiKj'^B,,o{p), X = 0,1, (21) 

where R^ is the Wigner rotation, L^, is the Lorentz boost from the baryon rest frame to a 
moving one, and the rotation matrix of the heavy quark spin D^^'^(R) in spinor representation 
is given by 

(J fjD'ci^RYj = S-\pQ,)S{A)S{p), (22) 

where 



siP) = \MA^^+ ""^ 



2m \ e{p) + mj 

is the usual Lorentz transformation matrix of the four-spinor. The rotation matrix D^{R) of 
the diquark with spin T is equal to D^{R}^) = 1 for the scalar diquark and D]^{R^) = R^ 
for the axial vector diquark. 



IV. FORM FACTORS AND SEMILEPTONIC DECAY RATES 

In this section we give the general parameterization of semileptonic decay matrix elements 
and the expressions for decay rates of heavy baryons with scalar and axial vector light 
diquarks. 

A. Heavy baryons with the scalar diquark 

The hadronic matrix elements for the semileptonic decay Aq -^ Ag/ are parameterized 
in terms of six invariant form factors: 



{Kq,{v\s')\V^\Aq{v,s)) = u^^Xv'.s) 
{Aq,{v\s')\A^\Kq{v,s)) = UA^,{v',s') 



UAq[V,S), 



Gi{w)^^ + G2{w)v^ + G^{w)v'^ 1^ukq{v, s), (23) 



where u\q{v, s) and u\ ,{y' ^ ^') ^^^ Dirac spinors of the initial and final baryon with four- 
velocities V and f', respectively; q = Ma v' — M^^v, and 

, Ml,+Ml -g^ 

W = V ■ V = ——z -^ . 

The helicity amplitudes are expressed in terms of these form factors '2l'] as 
1 



^)^ _ /oTJ" TF. /'„,, -I- 1 AT/' i\/r . _i_ i\/r. \ Tr^'^f„.,\ -L- i\/r /'„,, _i_ i ^ t^'^, 



^1/2 = -^j2MA^MAjwTmMA^±MA^,)J'^''^{w)±MAjw±l)J^, 



-V,A 



w) 



±Ma^(u;±1)^3^'"H], 
H^lt,i = -2^Ma^Ma^X^tI):F^'^{w), (24) 

where the upper(lower) sign corresponds to V{A) and T^ = Fi, Tf" = Gi {% = 1,2,3). 
Hy ^^ are the helicity amplitudes for weak transitions induced by vector {V) and axial 
vector (A) currents, where A' and Xw are the helicities of the final baryon and the virtual 
PF-boson, respectively. The amplitudes for negative values of the helicities can be obtained 
using the relation 

—X',—Xw ^' ^ ^w ' 

The total helicity amplitude for the V — A current is then given by 

TJ — IjV TjA 

J^\',Xw — ^X',Xw ~ ^X',Xw- 

The total differential decay rate 

^ = ^ + ^ (25) 

dw dw dw 

is expressed in terms of the partial rates for transversely (T) and longitudinally (L) polarized 
l^-bosons 



^ - (2;o^'^««'' — umT, — ^'^z'-^' +1^-1/2, -il J, 

where Gp is the Fermi coupling constant and Vqqi is the relevant CKM matrix element. 

The decay products in the semileptonic decay Aq — > Aq/(— > Avr) + W{-^ lu) are highly 
polarized. The polarization of the decay products is usually expressed through different 
asymmetry parameters l21] defined as follows: 

\U |2 I C7 |2 

-ni/2.1 — -n-l/2,-1 



a-L 



a 



I W 12 I 7T 12 

|-ni/2,o| — |-n-l/2,o| 

|ifi/2,oP + |^-iAoP' 

|-ni/2,l| — |-n-l/2,-l| 

\H^/2,i? + \H-y2,-i? + 2(|i/i/2,oP + |i^-i/2,on ' 



10 

// _ |-^l/2,l| — 1-^-1/2, -1 1 — 2(|iJ"i/2,o| + I -H"-^ 1/2,0 1 ) 

" ~ l^iAiP + \H-i/2,-i\^ + 2(|/7i/2,o|2 + l^-iAoH ' 

^ 2Re(iJ'_i/2, 0^1/2,1 + -^1/2,0^-1/2,-1) ^27^ 

'^ I f7 12 I I Z17 12 I \U 12 I I Z17 12' \ ') 

|-ni/2,ir + |-n-i/2,-ir + |-ni/2,or + |-n-i/2,or 

The average values of these asymmetry parameters ((ar), (ol), (a'), («"), (7)) are calculated 
by separately integrating the numerators and denominators in (j^7|) over w. The average 
longitudinal Ag/ polarization (Pl) can be expressed in terms of (ax), {(11) as 

^p^^^ (arHRM ^ fl=Ii. (28) 

B. Heavy baryons with the axial vector diquark 

The hadronic matrix elements for the semileptonic decay VLq -^ VLqi are parameterized 
in terms of six invariant form factors through expressions analogous to Eqs. (|23|). Then the 
helicity amplitudes and differential decay rates are given by Eqs. (jH)), (f^ with obvious 
mass replacements. 

The invariant parameterization for the semileptonic decay VLq -^ VL*q, reads: 

{n*Q,{v\s')\V^'\nQ{v,s)) = un'^x{v\s') [iViHA'^ + iV2Ht;V 

+K^{w)v\'^' + K^{w)g'^^']unQ{v,s), (29) 
where mq* ^^ is the Rarita-Schwinger spinor for the fig, which obeys 

l^un*^A'"^ *) = ^^qA^"^ ^)' v^'un^^^f.iv, s) = -f^un*^,f,{v, s) = 0. 
The helicity amplitudes are given by [2l| 

TiMn^ T Mo^,)(w ± m^'^'iw) + Mn^Jw^ - m^'^w) 

H^/t = \[l^/Mn^Mn^^,{^ T 1)[A/T(«;) - 2{w ± 1)<'^(«;)], 

H^jt^, = T^2Mn^MnliwTl)K'''M, (30) 

where again the upper (lower) sign corresponds to V{A) andAf^ = Ni, Aff^ = Ki {i = 1, 2, 3). 
The remaining helicity amplitudes can be obtained using the relation 

ttV,A _ _^ ttV,A 

Partial differential decay rates can be represented in the following form 



T ^P i9 O' i9 I i9 I i9 I i9i 

^ = (^I^QQ'I rn^^T [I^V2,ll +l^-l/2,-l| +|i^3/2,l| +|i^-3/2,-l|], 
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"^^ _ JfLF_|T/ |2__V \\TT iSiIrr |2] /oi \ 



V. SEMILEPTONIC DECAYS OF HEAVY BARYONS WITH THE SCALAR 

DIQUARK 

In the heavy quark hmit mq ^ oo {Q = b, c) the form factors (j^ can be expressed 
through the single Isgur-Wise function Cl""^) ll^ 

FiH = G^{w) = CH, 

F2{w) = Fsiw) = G2iw) = Gsiw) = 0. (32) 

At subleading order of the heavy quark expansion two additional types of contributions arise 
221 . The first one parameterizes I/itiq corrections to the HQET current and is proportional 



to the product of the parameter A = Maq — mq, which is the difference of the baryon and 
heavy quark masses in the infinitely heavy quark limit, and the leading order Isgur-Wise 
function C('^)- The second one comes from the kinetic energy term in l/mq correction to 
the HQET Lagrangian and introduces the additional function xi'^)- Therefore the form 
factors are given by }22|| 



F^{w) = Ciw) + ( 7T^ + ^] [2x{w) + Ciw)] , 



v2mQ 
G,{w) = C{w)+(-^+ ^ 



2mQ 2mQi 



If — 1 
w + 1 



F,(y,) = -G3W = -5^^C(^). (33) 

To calculate these semileptonic decay form factors in our model we substitute the vertex 
functions F^^^ (jl8|) and F*^^) (|19p in the weak current matrix element (J17p between Aq and Ag/ 
baryons. It is important to take into account the relativistic transformation of the baryon 
wave functions (PT|) in this matrix element. The resulting structure of the decay matrix 
element is rather complicated, because it is necessary to integrate both over d^p and d^q. The 
5 function in expression (fTH|) for F*^^^ permits us to perform one of these integrations and thus 
this contribution can be easily calculated. The calculation of the contribution of the vertex 
function F*^^^ (fT^ is more difficult, since here, instead of a 5 function, we have a complicated 
structure, containing the heavy-quark-light-diquark interaction potential. Nevertheless, the 
application of the heavy quark l/mg expansion considerably simplifies the calculation. We 
carry out this expansion up to the first order. Then we use the quasipotential equation 
to perform one of the integrations in the decay matrix element. The vertex function F*^^) 
provides the leading order contribution, while F*^^-* contributes already at the subleading 
order. The resulting expressions for the semileptonic decay Ag -^ Ag/ form factors up to 
subleading order in l/rrig are then given by 
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A 



2mo' w — 1 2m. 



GiH = CH + 




iw+l) 
Q 

w + I 



X[w) 



-4{1 - e){l + k)- wx{w), 

2mn 



Fo(w) 



A 



C(^) 



G2{W) 

Fs{w) 



2mQi w + 1 
A 2 



A 



1 A 

+ ::: W 



2mn' w — 1 2m 



2m,Qr w + 1 
-Gsiw) = 



CM -4(1 -£)(! + «:) 



A 



Xiw) 
1 



2m.ni w — \ 



X{w) 



2mo w + 1 2mo 



(34) 



where the leading order Isgur-Wise function of heavy baryons 



CH 



hm 



(2vr) 



^Ao,| P + 2ed(p)W^^ eA I ^Ag(p), 



3 '^Q 



W +1 



(35) 



and the subleading function 



x{w) 



w — 1 



hm 



W + 1 rriQ^oo J (27r 



'''^*.ip + 2.,(p)J^ e.) ^^?M*,^(p), (36) 



\3 '^Q' 



W+l 



2A 



here ga 



A/v A^ is the unit vector in the direction of A = Ma v' — Ma^v. It is 
important to note that in our model the expressions for the Isgur-Wise functions C{w) (j35p 
and x{w) (|36|) are determined in the whole kinematic range accessible in the semileptonic 
decays in terms of the overlap integrals of the heavy baryon wave functions, which are 
known from the baryon mass spectrum calculations. Therefore we do not need to make any 
assumptions about the baryon wave functions or/and extrapolate our form factors from the 
single kinematic point, as it was done in most of previous calculations. 

For {l — e){l + K) = the HQET results (j^Hj) are reproduced. This can be achieved either 
setting e = 1 (pure scalar confinement) or k = — 1. In our model we need a vector confining 
contribution (see Sec. IHl) and therefore use the latter option. This consideration gives us an 
additional justification, based on the HQET, for fixing one of the main parameters of the 
model K."^ In the heavy quark limit the wave functions of the initial \1/aq and final baryon 
\1/a„, coincide, and thus the HQET normalization condition ((1) = 1 is exactly reproduced. 
The subleading function x(^) vanishes for w = 1. These functions, calculated with model 
wave functions for A^ and A^ baryons, are plotted in Figs. El IH The function xi'^) is very 
small in the whole accessible kinematic range, since it is roughly proportional to the ratio 
of the heavy baryon binding energy to the baryon mass. 



^ It is important to note that the same value of k is needed_to get agreement with the HQET structure of 
the first order l/mq corrections for B -^ D'^*^ev decays 
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FIG. 3: The Isgur-Wise function C(^) foi^ the Af, — > Acez^ semileptonic decay. 



TABLE III: 


Parameters of the Isgur 


-Wise functions 


for the Afo - 


-^A, 


ei^ 


and Hfo - 


^ Hcez^ decays. 


Decay 


A (GeV) 


2 

Pc 


CC 






P^. 


Cx 


Afe ^ Accz^ 
Eb -^ Eceu 


0.764 
0.970 


1.70 

2.27 


2.39 

3.87 






0.053 
0.045 


0.029 
0.036 



Near the zero recoil point of the final baryon w = 1 these functions can be approximated 



by 






1-pI{w 

2 



l) + cc(^-l)' + - 



pJw-l) + Cx{w 



l? + 



(37) 



where p^ = —[d({w) / dw]ii;=i is the slope and 2c^ = [(f({w)/(fw]w=i is the curvature of 
the Isgur-Wise functions, which are given in Table ITTTl The values of ({1) for transitions 
between physical Af, (S^) and A^ (S^) baryons are slightly different (by ~ 0.5%) from the 
heavy quark limit value 1 due to the distinction of the A;, (S;,) and Ac (Sc) baryon wave 
functions, calculated for finite values of the heavy quark masses. 

Our model predictions for the form factors Fi{w) and Gi{w) {i = 1, 2, 3) for the Af, -^ AcCu 
semileptonic decay are plotted in Fig. El The corresponding differential decay distributions 
calculated both with inclusion of first order heavy quark corrections and in the heavy quark 
limit are plotted in Fig. [Ul 

The Aft -^ Ac differential decay rate near zero recoil I2i 



1 dT{Ab 
lim . I 



Accz/) 



dw 



Gp\Vcb\ 

47]-3 



-M1JMa,-MaJ2|Gi(1)P 



(38) 
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FIG. 4: The subleading function x(^) for the Af, -^ K^ev semileptonic decay. 



is governed by the square of the axial current form factor Gi, which near this point has the 
following expansion 

Gi{w) = 1 - p\w - 1) + c(«; - 1)2 + ■ ■ ■ , (39) 

where in our model with the inclusion of the first order heavy quark corrections (|34j) 



p^ = 1.51, 



and c = 2.03. 



This value of the slope parameter of the Af,-baryon decay form factor is in agreement with 
the recent experimental value obtained by the DELPHI Collaboration J23| 



p2 = 2.03 ± 0.46+? 



00 



and lattice QCD J2J] estimate 

p2 = 1.1 ±1.0. 

Our prediction for the branching ratio of the semileptonic decay A^ -^ A^eu for |V"cfe| 
0.041 and ta, = 1.23 x lO^^^g jj^ 

Br'^'°\Ab^AJu) =6.9% 



is in agreement with available experimental data 



5-0tlit\f) % DELPHI [23] 

8.1±1.2li;^±4.3)% CDF [25] 



and the PDG branching ratio |1^ 

5r'^P(Afe -^ AJu + anything) = (9.1 ± 2.1)%. 



(40) 



(41) 
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FIG. 5: Semileptonic decay form factors for A^ -^ AcCu. 

TABLE IV: Theoretical predictions for semileptonic decay rates T (in lO^'^s"^) and averaged asym- 
metries for Ah — > AcCv and H^ — > H^ei^ for \Vcb\ = 0.041. Branching ratios Br (in %) are calculated 



using experimental mean 


values 


u 


:or the life times 


TA, = 1.23 


X 10-12s 


and thj, 


= 1.39x 


10-i2s. 


Decay T Br 


Tl 


Tt 


R 


(ot) 


(«l) 


{Pl) 


(«') 


{a") 


(7) 


in mq -^ oo limit 




















Afc -^ Accu 5.02 6.2 


3.08 


1.94 


1.59 


-0.483 


-0.928 


-0.756 


-0.116 


-0.521 


0.562 


Eb -^ Eccu 4.64 6.4 


2.79 


1.85 


1.51 


-0.455 


-0.920 


-0.735 


-0.113 


-0.503 


0.587 


with l/mg corrections 




















Afo -^ Acev 5.64 6.9 


3.48 


2.16 


1.61 


-0.600 


-0.940 


-0.810 


-0.142 


-0.527 


0.494 


Hb ^ HcCi^ 5.29 7.4 


3.21 


2.08 


1.54 


-0.597 


-0.935 


-0.802 


-0.146 


-0.510 


0.505 



Predictions of our model for the semileptonic decay rates (j^UI) and averaged asymmetries 
()27|) and ()28p for Aj, -^ AcCz/ and S^ — > ScCz/ decays, both in the heavy quark limit and with 
inclusion of first order l/niQ corrections, are given in Table HVl In decay rate calculations we 
used for the Sh mass the value from Table ITTl and for other masses their experimental values 
[l9|. Comparing results for the decay rates with and without first order I/itiq corrections 
we see that the inclusion of the subleading terms leads to a relatively small ~ 14% increase 
of the total decay rates. Therefore, one can expect that higher order corrections should be 
small, and thus their account cannot substantially change the leading order predictions. 
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FIG. 6: Differential decay rates dT/dw for the A^ -^ K^ev semileptonic decay. Solid lines show 
decay rates including first order l/mq corrections. Dashed lines correspond to decay rates in the 
heavy quark limit. 



VI. SEMILEPTONIC DECAYS OF HEAVY BARYONS WITH THE AXIAL 

VECTOR DIQUARK 

In the heavy quark limit uiq -^ oo the decay matrix element ()29p is reduced to 12, ll5| 
(fi«K, s')\hi?'^Th^^^\nQ{v, s)) = 55>', s')TB^^{v, s)[-g^''Ciiw) + t;V'^C2(^)], (42) 



where 



B^^iv,s) = -^ 



{l^^ + v^,)-f5UnQ{v,s), 



Bf,'^{v,s) =un* ^^{v,s). 



(43) 



The structure of the leading order in l/mq corrections, which in the HQET can be param- 
eterized in terms of five additional functions, can be found in Ref . [l^ . 

In our model the corresponding semileptonic decay matrix element can be calculated using 
the same procedure as in the previous section. However such calculation is considerably more 
cumbersome (especially for the r'-^-' contribution), since now the spectator light diquark has 
spin equal to 1. Taking into account that in the Aq baryon decays contributions of l/mq 
corrections are rather small, we expect that in the case of Qq baryon decays such corrections 
should be also relatively small. At present no bottom baryons with axial vector diquark have 
been observed yet, and, when observed, their semileptonic decays will be difficult to measure. 
Therefore it seems reasonable to limit our analysis here to the leading order of the heavy 
quark expansion. In the heavy quark limit only the lowest order vertex function F*^^^ (I18p 
contributes to the decay matrix element ()17|) . The resulting expressions for the weak decay 
matrix element exactly satisfy the HQET relation (jl^ and allow us to determine the Isgur- 
Wise functions Ci('^) ^^^ C2{w) in the whole accessible kinematic range through the overlap 
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FIG. 7: The Isgur-Wise functions Ci(^) and C2{w) for the Jl^ -^ Qh ev semileptonic decay. 



W. 



integrals of the baryon wave functions. They are given by 



C2(^) 



hm 



(2vr)^ 



r^n«, P + 2ed(p)i 



If 



w + 1 



eA^f^cslP) 



QVi'y) 



W + l 



Cl(«^), 



(44) 
(45) 



where e^ = A/vA^ is the unit vector in the direction of A = M^ v' — Mq^v. The 
relation ()45|) follows from the relativistic spin transformation ()21|) of the spectator axial 
vector diquark. A similar relation was obtained also in Ref. |3]. The Isgur-Wise functions 
are plotted in Fig. [7| 

Near the zero recoil point w = 1 the Isgur-Wise functions can again be approximated by 



Uw) = 0(1) - pUw - 1) + c^w - 1)2 + 



(46) 



where Ci(l) = 1 and C2(l) = 1/2; Pq = —[d(i{w)/dw]w=i is the slope and 2cq = 
[(P(i/d'^w]w=i is the curvature of the Isgur-Wise functions, which are given in Table M 

The invariant form factors in the heavy quark limit can be expressed, using relation (j45|) . 
in terms of the Isgur-Wise function Ci('"^) as follows 

F2{w) = l-^Ci{w), 



G2{w) 

Ni{w) 



3w + l 
Gs{w) = 0; 

-N,{w) 



K.iw) 



\/3w + l 



(iM 
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TABLE V: Parameters of the Isgur-Wise functions for the Sj 






•J c e.v and fi^ 



f^c ei^ decays. 












Decay 


A (GeV) 


p?, 


cci 


p?. 


CC2 


Eft ^ S^*^ezy 


0.942 


2.17 


3.62 


1.34 


2.44 




1.082 


2.61 


4.93 


1.55 


3.19 


rjfe ^ f7i*^ez^ 


1.208 


2.99 


6.21 


1.74 


3.91 



ro 



^ 2 




1.1 



1.2 



1.3 



tf; 



FIG. 8: Differential decay rates dV /dw for the Of, ^ Ocei/ semileptonic decay. 



N2(w) = Ki(w 



K2(w) = 0. 



(47) 



The differential decay rates dV /dw for ^2^ -^ Vt'^'eu semileptonic decays are plotted in 
Figs. IHl and El The decay rates of bottom baryons with the axial vector diquark, calculated 
in the heavy quark limit using expressions (|!?T|l , are given in Table IVII For masses of bottom 
baryons and the f2* we used the values from Table |n] and for other charmed baryons we used 
experimental mass values jl9|. 



VII. DISCUSSION AND CONCLUSIONS 



The cornparison of our model predictions with other theoretical calculations 



i ne cornparison oi our moaei preaictions witn otner tneoreticai calculations }ci |zhi Mi , 
l28l . I29L I.SflL l3ll . l32| | is given in Table IVIII In nonrelativistic quark models J3, 120, |22| form 
factors of the heavy baryon decays are evaluated at the single kinematic point of zero re- 
coil and then different form factor parameterizations (pole, dipole) are used for decay rate 
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FIG. 9: Differential decay rates dV/dw for the $7^ — > f]*ez^ semileptonic decay. 



TABLE VI: Theoretical predictions for semileptonic decay rates T (in 10^'^s 
with the axial vector diquark in the heavy quark limit for \Vcb\ = 0.041. 



of bottom baryons 



Decay 


r 


Tl 


Tr 


R 


Sfi -^ TicCv 


1.44 


1.23 


0.21 


5.87 


E', -. E'^eu 


1.34 


1.14 


0.20 


5.59 


Qb -^ ^cev 


1.29 


1.09 


0.20 


5.31 


Sfc ^ T^leu 


3.23 


1.61 


1.62 


0.99 


Ei ^ E*eu 


3.09 


1.52 


1.57 


0.97 


Qb —>■ ^leu 


3.03 


1.48 


1.55 


0.95 



calculations. The relativistic three-quark model [28[, Bethe-Salpeter model j29|] and light- 
front constituent quark model |30| assume Gaussian wave functions for heavy baryons. The 
authors of the recent nonrelativistic quark model fSl^ use for the form factor evaluations 
the set of variational wave functions, obtained from baryon spectra calculations without 
employing the quark-diquark approximation. Finally, Ref. |32] presents the recent QCD 
sum rule prediction. Calculations of Refs. J3|, |2^ |29| are done in the heavy quark limit only, 
while the rest include first order l/mq corrections for the decays of A-type baryons. From 
Table IVTTl we see that all theoretical models give close predictions for the semileptonic decays 
of heavy baryons with scalar diquark (A;, -^ AcCu and S;, -^ Eceu), which are consistent 
with the available experimental data ()40j) and ()41|) for the A;, — > AcCz/ semileptonic decay. 
The results for averaged asymmetries of these decays (see Table IIV|) are also close in most 
of the considered approaches. Thus one can conclude that the precise measurement of the 
semileptonic A^ -^ A^eu decay rate will allow an accurate determination of the CKM matrix 
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TABLE VII: Comparison of different theoretical predictions for semileptonic decay rates T (in 
lO^'^s"^) of bottom baryons. 



Decay 


this work 


[26j 


[27] 


[3] 


[28] 


[29] 


[30] 


[31] 


[32] 


A(, -^ KcCv 


5.64 


5.9 


5.1 


5.14 


5.39 


6.09 


5.08 ±1.3 


5.82 


5.4 ±0.4 


Eb -^ Eccu 


5.29 


7.2 


5.3 


5.21 


5.27 


6.42 


5.68 ±1.5 


4.98 




T,b -^ TicSv 


1.44 


4.3 






2.23 


1.65 








Ei ^ E'^eu 


1.34 


















^b -^ ^c^y 


1.29 


5.4 


2.3 


1.52 


1.87 


1.81 








T.b -^ T^lev 


3.23 








4.56 


3.75 








^6 -^ "^1^^ 


3.09 


















nb -^ nieu 


3.03 






3.41 


4.01 


4.13 









element Vcb with small theoretical uncertainties. 

All predictions for heavy baryon decays with the axial vector diquark listed in Table IVllI 
were obtained in the heavy quark limit. Here the differences between predictions are larger. 
The nonrelativistic quark model J2h| gives for these decay rates values more than two times 
larger than other estimates. Our model values for these decay rates are the lowest ones. 
Among the relativistic quark models the closest to our predictions is given in [29[. Unfor- 
tunately, it will be difficult to measure such decays experimentally. Only Qb (which has 
not been observed yet) will decay predominantly weakly and thus has sizable semileptonic 
branching fractions, since a scalar ss diquark is forbidden by the Pauli principle. All other 
baryons with the axial vector diquark will decay predominantly strongly or electromagneti- 
cally and thus their weak branching ratios will be very small. 

In summary, in this paper we calculated the semileptonic decay rates of heavy baryons in 
the framework of the relativistic quark model. Heavy baryons were considered in the heavy- 
quark-light-diquark approximation. The baryon wave functions were obtained previously in 
the process of the heavy baryon mass spectrum calculations. In our approach the spectator 
diquark is not treated as a point-like object. The relatively large diquark size is taken 
into account by calculating the diquark-gluon form factor as the overlap integral of the 
diquark wave functions. The matrix element of the weak current between baryon states was 
considered using the quasipotential approach. The relativistic transformation of the baryon 
wave functions from the rest reference frame to the moving one as well as the negative energy 
contributions to the decay matrix elements were explicitly taken into account. To simplify 
calculations and in order to compare with model-independent predictions of HQET the heavy 
quark expansion was applied up to subleading order for heavy baryon decays with a scalar 
light diquark. It was shown that all HQET relations in the leading and subleading order 
are exactly satisfied in our model if the long-range chromomagnetic interaction vanishes 
(k = —1) in accord with our previous analysis of heavy meson decays. The leading and 
subleading Isgur-Wise functions were determined in terms of the overlap integrals of baryon 
wave functions. It was found by explicit calculation that the additional subleading function 
x{w), arising from the kinetic energy term in the HQET Lagrangian, is negligibly small in 
the whole kinematic range. Decay rates as well as different averaged asymmetries both with 
and without l/mg corrections were calculated. Moreover, it was shown that the subleading 
terms in the heavy quark expansion modify the results for decay rates by ~ 14%. Thus one 



21 

can expect that the influence of higher order corrections should be small. 

The decays of heavy baryons with the axial vector diquark were considered in the heavy 
quark limit. All HQET relations are exactly satisfied in our model. The corresponding 
Isgur-Wise functions were determined in terms of the overlap integrals of the baryon wave 
functions. It was found that the relativistic transformation of the axial vector diquark spin 
leads to the relation (jlKjl between baryon Isgur-Wise functions Ci(''^) ^^^ C2{w). 

The calculated decay rates of heavy baryons were compared with the results of other 
theoretical approaches and available experimental data. One of the main advantages of 
our model is that it allows one to calculate consistently the heavy baryon wave functions 
from the consideration of the spectroscopy and then determine through the wave function 
overlap integrals the baryonic Isgur-Wise functions in the whole kinematic range accessible 
in semileptonic decays. Thus we do not need to make any assumptions about the form of 
the baryon wave functions or/and extrapolate the form factors from one point to the whole 
kinematic region using some ad hoc ansatz. No additional free parameters were introduced 
in our calculations. As it was pointed out above, we also consistently include relativistic 
effects. All this makes the presented results sufficiently accurate and reliable. 
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